Evidence for tt + tt scattering in p + p collisions at 
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Abstract. A p(770)° mass shift of about -40 MeV/c 2 was measured in p+p collisions 
a t -J s nn = 200 GeV at RHIC. Previous mass shifts have been observed at CERN- 
LEBC-EHS and CERN-LEP. We will show that phase space does not account for 
the /9 (770) mass shift measured at RHIC, CERN-LEBC-EHS and CERN-LEP and 
conclude that there are significant scattering interactions in p + p collisions. 
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1. Introduction 

The p° was measured via its hadronic decay channel in minimum bias p + p and Au+Au 
collisions at RHIC and a mass shift of -40 and -70 MeV/c 2 of the position of the p° 
was observed, respectively. While no explanations were explicitly attributed to the 
mass shift in p + p collisions, the possible explanations for the apparent modification of 
the p° meson properties in Au+Au collisions were attributed to dynamical interactions 
with the surrounding matter, interference between various 7r + 7r~ scattering channels, 
phase space distortions due to the rescattering of pions forming p° and Bose-Einstein 
correlations between p° decay daughters and pions in the surrounding matter [T|. 

The p° meson measured in the dilepton channel probes all stages of the system 
formed in relativistic heavy-ion collisions because the dileptons have negligible final 
state interactions with the hadronic environment. Heavy-ion experiments at CERN 
show an enhanced dilepton production cross section in the invariant mass range of 200- 
600 MeV/c 2 , showing that the p° is broadened rather than shifted [21 E]- The hadronic 
decay measurement at RHIC |T], p(770)° — > 7r + 7r+ was the first of its kind in heavy-ion 
collisions. Since the p° lifetime of cr = 1.3 fm is small with respect to the lifetime of 
the system formed in Au+Au collisions, the p° meson is expected to decay, regenerate, 
and rescatter all the way through kinetic freeze-out. Therefore, the measured p° mass 
at RHIC should reflect conditions at the late stages of the collisions E] ■ However, its 
has been shown that the p° width goes to zero with the p° mass dropping near the chiral 
phase transition, which is called the the vector manifestation and it can occur near T c 
in the hot and dense matter In this scenario, the p° lifetime becomes very large and 
the p° produced at the early stages can also be measured in the hadronic channel. 

The modification of the p° properties in heavy-ion collisions has been expected 
0], contrary to the modifications of the p° properties in p + p collisions. Similar 
modifications of that measured in p + p collisions at RHIC has been observed before 
at CERN-LEBC-EHS and CERN-LEP. Until now, the mass shift measured in p + p at 
CERN-LEBC-EHS has been attributed to phase space. The process of hadron-hadron 
collisions at low pr is not well understood. These collisions are viewed as a complicated 
interaction, where a perturbative treatment of the interaction dynamics is not possible. 
As a consequence, only phenomenological models with assumptions about the partonic 
subprocess dynamics can be used to interpret the data. By studying soft hadronic 
interactions we may be able to explain the interaction dynamics of hadron collisions 
and possibly understand the parton hadronization mechanism, which is important in 
jet studies. One way of probing the dynamics of hadron-hadron collisions is to study 
the production of resonances, in particular, the p° meson. 

In this paper, we will show that phase space does not account for the p(770)° mass 
shift measured at RHIC 0, CERN-LEBC-EHS [B] and CERN-LEP QUJ El CU and 
conclude that there are significant scattering interactions in p + p collisions. 
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2. Discussion p mass average from the Particle Data Group (PDG) 

We will first discuss the p mass average from the PDG |14j . The p° mass average 775.8 
± 0.5 MeV/c 2 from e + e~ was obtained from either e + e~ — > 7r + 7r~ or e + e~ — > 7r + 7r~7r°. 
This means that the p° mass average was obtained from exclusive leptonic reactions. 
Similarly, the p ± mass average 775.5 ± 0.5 MeV/c 2 was also was obtained from exclusive 
leptonic reactions. 

The p averages reported by the PDG from reactions other than leptonic interactions 
are systematic lower than the value obtained from leptonic exclusive interactions by ~ 10 
MeV/c 2 [Hj. The p production in these hadronic reactions are inclusive and exclusive. 
In the case of inclusive productions, the phase space was take into account when the p 
mass was measured (e.g. jH]). 

These observations lead us to conclude that the p mass depends on specific 
interactions, e.g. whether the p is produced in inclusive or exclusive reactions. Since a 
leptonic reaction and exclusive measurement of the p° lead to a negligible modification 
of any kind of the p° mass, the average 775.8 ± 0.5 MeV/c 2 from e + e~ should 
correspond to the p° mass in the vacuum. 



3. p° mass shifts at RHIC, CERN-LEBC-EHS, and CERN-LEP 

The p° was measured in the hadronic decay channel p° — > tx + tx~ at RHIC, CERN- 
LEBC-EHS, and CERN-LEP in inclusive production. At RHIC, the STAR collaboration 
measured the p° at ^/s NN = 200 GeV at midrapidity (\y\ < 0.5) and observed mass shifts 
of the position of the p° peak of about -40 MeV/c 2 and -70 MeV/c 2 in minimum bias p+p 
and peripheral Au+Au collisions, respectively p. The invariant mass distributions from 
PP is shown in Fig. ^ The solid black line in Fig. Q is the sum of all the contributions 
in the hadronic cocktail. The K® was fit to a Gaussian (dotted line). The u (light grey 
line) and K*(892)° (dash-dotted line) shapes were obtained from the HIJING event 
generator [T^], with the kaon being misidentified as a pion in the case of the K*°. The 
p°(770) (dashed line), the / (980) (dotted line) and the / 2 (1270) (dark grey line) were fit 
by relativistic Breit-Wigner functions (BW) times the Boltzmann factor (PS) to account 
for phase space |3 HBl CZ1 HH1 CHI 12H 122], 

BW= M - MoF (1) 

(M 2 - M 2 J 2 + M 2 r 2 1 } 

r _r nX ii x ^l-Aml (2m)/2 

r " r ° X M^ X[ M 2 -4m 2j (2) 



oc_ xex -V^+Pt 



where Mq and To are the natural resonance mass and width, respectively PUE]. The 
masses of Kg, p°, f , and / 2 were free parameters in the fit, and the widths of p°, f and 
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Figure 1. The raw ir + ir~ invariant mass distributions after subtraction of the like-sign 
reference distribution for minimum bias p + p (top) and peripheral Au+Au (bottom) 
interactions measured by STAR. For details see 



f 2 were fixed according to [14 . The uncorrected yields of K s , p°, u>, fo, and f 2 were free 
parameters in the fit while the K*° fraction was fixed according to the K*(892) ^ttK 
measurement (22]. The p°, u, K*°, f , and f 2 distributions were corrected for the 
detector acceptance and efficiency determined from a detailed simulation of the TPC 
response using GEANT pQ. The number of degrees of freedom (d.o.f.) from the fits 
was 196 and the typical xVd.o.f. was 1.4. The p° mass obtained from the BWxPS fit 
is depicted in Fig. El where it is clear that the phase space does not account for the 
measured mass shifts of the position of the p° peak. 

At CERN-LEBC-EHS, NA27 measured the p° in minimum bias p + p at = 
27.5 GeV for xp > 0, where xp is the ratio between the longitudinal momentum and 
the maximum momentum of the meson, and reported a mass of 762.6 ± 2.6 MeV/c 2 
jH]. The p° signal can be seen even before background subtraction. The invariant tt + tt~ 
mass distribution was fit to the BW(NA27)xPS(NA27) plus a background function [B] 

BW(NA27) = JM^m X 

(4) 



{Ml - M^y + M 2 P 



r nX m q :;[ M i- 4m ' i(2W)/ 2 (5) 
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Figure 2. The p° mass as a function of for minimum bias p+p (filled circles), high 
multiplicity p+p (open triangles), and peripheral Au+Au (filled squares) collisions 
measured by STAR. The error bars indicate the systematic uncertainty. Statistical 
errors are negligible. The p° mass was obtained by fitting the data to the BWxPS 
functional form described in pQ . The dashed lines represent the average of the p° mass 
measured in e + e~ 14_. The shaded areas indicate the p° mass measured in p + p 
collisions [H]- The open triangles have been shifted downward on the abscissa by 50 
McV / c for clarity. 



PS(NA 2 T, = BG = ( Z)^uT 1 ^ * 

exp(-a3[%M]^- ai[ %M]) (6) 

where ai, 02,03, and 04 are free parameters in the fit. In this analysis, the phase 
space function (PS(NA27)) used is the same as the combinatorial background (BG). 
The invariant tt + tt~ mass distribution after subtraction of the mixed-event reference 
distribution is shown in Fig. El The vertical dashed line represent the average of the 
p° mass 775.8 ± 0.5 MeV/c 2 measured in e + e~ fll]. The vertical solid line is the p° 
mass 762.6 ± 2.6 MeV/c 2 reported by NA27 |8 r As shown in Fig. El the position of 
the p° peak is shifted by ~ 30 MeV/c 2 compared to the p° mass in the vacuum 775.8 ± 
0.5 MeV/c 2 [Hj. The shift of the p° peak can be quantified by fitting the NA27 tt+tt" 
mass distribution after subtraction of the mixed-event reference distribution to the BW 
function (equation^). The fit is depicted in Fig. HJand the value of the maximum of the 
distribution obtained from the fit is 747.6 ± 2.0 MeV/c 2 with x 2 /ndf = 1.9. If equation 
0]is used instead, the value of the maximum of the distribution is 754.3 ±2.1 MeV/c 2 
with x 2 / n df = 1-2. The difference of a few MeV between the two values is due to the 
different BW functions used in the analyzes. 

As mentioned, NA27 obtained the p° mass by fitting the invariant 7r + 7r - mass 
distribution to the BWxPS function, and they reported a mass of 762.6 ± 2.6 MeV/c 2 , 
which is ~ 10 MeV/c 2 lower than the p° mass in the vacuum. Ideally, the PS factor 
should have accounted for the shift on the p° peak, and the mass obtained from the fit 
should have agreed with the p° mass in the vacuum. However, just like in the STAR 
measurement, this was not the case, since the phase space did not account for the mass 
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Figure 3. The invariant ir + ir~ mass distribution after background subtraction 
(equation for minimum bias p + p collisions measured by NA27. For details see 
[Hj. The vertical dashed line represent the average of the p° mass 775.8 ± 0.5 MeV/c 2 
measured in e + e~ 14 . The vertical solid line is the p° mass 762.6 ± 2.6 MeV/c 2 
reported by NA27 0. 
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Figure 4. The invariant ir + ir~ mass distribution after background subtraction 
(equation |SJ) for minimum bias p + p collisions measured by NA27 |Hj fit to the BW 
function (equation 0. For details see |S]. 



shift on the position of the p° peak. 

At CERN-LEP, OPAL, ALEPH and DELPHI measured the p° in inclusive e + e~ 
reactions at yfs = 90 GeV Ell HU E2] ■ Even though OPAL never reported the value 
of the p° mass, OPAL reported a shift on the position of the p° peak by ~ 70 MeV/c 2 
at low x p , where x p is the ratio between the meson and the beam energies, and no shift 
at high x p (x p ~ 1) jniHD]- OPAL also reported a shift in the position of the peak 
from -10 to -30 MeV/c 2 , which was consistent with the p° measurement ^T]. ALEPH 
reported the same shift on the position of p° peak observed by OPAL • Both OPAL 
and ALEPH used the like-sign method to subtract the background. DELPHI fit the 
raw invariant 7r + 7r~ mass distribution to the (BWxPS + BG) function 

PS(DELPHI) = BG = (M — M th )~ n x 
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Figure 5. The pr distribution at \y \ < 0.5 for minimum bias p + p collisions |J. The 
black circles are the data, the dashed line is the power-law fit and the solid line is the 
exponential fit. See pQ for the fit functions. The errors shown are statistical only and 
smaller than the symbols. 



exp( 72 M + 73 M 2 + 74 M 3 + 75M 4 ), (7) 

where Mfh is the threshold mass and BW is a relativistic Breit-Wigner function, for x p > 
0.05 and reported a p° mass of 757 ± 2 MeV/c 2 [T3] . which is five standard deviations 
below the p° mass in the vacuum (775.8 ± 0.5 MeV/c 2 ). As one can see, similarly 
to NA27, DELPHI assumed that the phase space was described by the background 
function. Bose-Einstein correlations were used to describe the shift on the position of 
p° peak. However; high (even unphysical) chaocity parameters (A ~ 2.5) were needed 

p cni nil 

4. Phase space in p + p collisions 

In p+p collisions, most models assume that particles are born at hadronization according 
to phase space without any final state interaction. In multiparticle production processes, 
single-inclusive (e.g. p + p), invariant particle spectra are typically exponential in p^ 
[21]. The exponential behavior does not require final state interactions and it can be 
due to phase space population at hadronization. The slope parameter in p + p collisions 
are independent of the particle species [2H]. At RHIC, the p° spectra in minimum bias 
and high multiplicity p+p are exponential in p? up to 1.1 GeV/c 2 with slope parameters 
of ~ 180 MeV. The p° spectrum in minimum bias p + p pQ is depicted in Fig. 03 For 
reference, the slope parameter of tt~ is ~ 160 MeV [2E]- Note that the slope parameter 
in p + p is independent of the particle rest mass and shows m T scaling. These results 
are also independent of the multiplicity. 

However, the p° mass measured in p + p at RHIC is multiplicity dependent (the 
mass shift in high multiplicity is higher than in minimum bias p + p collisions) pQ, 
which is opposite to the slope parameter that is multiplicity independent. Note that the 
systematic errors are correlate among themselves and between the two measurements. 
We will demonstrate that the mass should be independent of px as well as multiplicity 
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at hadronization without final state interactions. According to quantum mechanics, a 
resonance at rest is described by the wave function, 

t) oc exp( l ^ ot ) x exp(— ^) (8) 

where r is the lifetime and Eq is the energy at rest. The probability amplitude of the 
resonance decay can be written as 

roo —iEt 

V(E,x)<x (eM-~—))*V(x,t)dt (9) 
Jo n 

^(E,x)cx 71E1 ^ r *(x) (10) 

P(E)k\ ^ e _1_ i | a / \9(x)\*dx (11 
P(£) oc V J1 (12) 

where T = fi/r, and the probability amplitude is a non-relativistic Breit-Wigner 
distribution. When the energy conservation law is imposed in the partition of string 
fragmentation into multiple hadrons, a phase space factor, similar to the Boltzmann 
factor in thermal distribution, has to be added to equation Replacing E by the 
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invariant mass (M), the phase space is exp(^p 1 ), where uit equals \J M 2 + p\ and T 
equals 160 MeV. Equation^] is then rewritten as 



P(M,p T )oc ^ ^xexp(^^). (13) 

1 ,V ' (M-M ) 2 + ^ T y 1 

However; as discussed previously, the phase space does not describe the mass shift of 
the p meson measured at RHIC, CERN, and LEP. 

Most event generators (e.g. PYTHIA [27] and HIJING |28j ) create resonances 
according to a non-relativistic Breit-Wigner function at a given p? (equation H!2|) . 
However, we just showed that a more reasonable way to produce resonances is to use 
equation [T31 

Since the phase space of a non interacting multiparticle state cannot explain the 
distortion of the p° line shape, we can conclude that the phase space in p + p collisions 
also accounts for hadrons scattering and forming resonances. In the case of the p°, 
7r + 7r~ — » p° — > 7T + 7r _ . This can be pictured in the string fragmentation particle 
production scenario, where the string breaks several times, two pions are formed, they 
scatter, and form a p . Such interactions are significant and modify the p° spectral 
shape in e + e~ and p + p from a relativistic p-wave Breit-Wigner function. 

Using the p + p data from the top plot of Fig. and the rescattering formalism of 
[2H], we can test the ideas of a mass shift having a 7r + 7r~ scattering component. The 
di-pion production is given by equation 21 of [2H]- For the p° which is p-wave I = 1 
equation 21 becomes 

|T| 2 = |£| 2 ^^ + ^ letting) + PS'«w(*i)| 2 (14) 
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Figure 6. Fit to the data using equation 1141 fsolid line), the direct p° contribution 
term times the phase space (dashed line), the rescattering term (dash-dotted line) and 
the sum of the both (dashed line). 



The 5i is the p-wave phase shift (p°) and sm p ^^ in our case becomes the Breit-Wigner 
times phase space equation[T31 \D\ 2 becomes the direct production of p° for the pr range 
600 to 800 MeV/c. \A\ 2 is the phase space overlap of di-pions in the £ = 1 partial waves. 
The pions emerge from a close encounter in a defined quantum state with a random 
phase. The emerging pions can re-interact or re-scatter through the p-wave quantum 
state or a phase shift. The phase space overlap comes from sampling the 7r spectrum 
from p + p collisions, where the sum of the 7r + 7r~ has the correct Pt range. The di- 
pion mass spectrum for this j>r range falls with increasing di-pion invariant mass. The 
variable a is related to the real part of the di-pion rescattering diagram and measures 
the range of interaction [29J 

2 

a =(1.0 -^) (15) 



r, 



The r is the radius of rescattering in fm and r is 1.0 fm or the limiting range of the 
strong interaction. The PS' is the phase space factor 

where q is the momentum of the pions in the center-of-mass and go is the size of 
interaction for the p° system and we use 1.0 fm distance to set go at 200 MeV/c. 

In order the fit the p + p data we need to add the rest of the cocktail of resonances 
used by STAR [T|. Such a fit is shown in Fig. |U] (solid line), where the direct p° term 
time the phase space (dotted line) and the rescattered term (dash-dotted line) plus the 
sum (dashed line) are also depicted . In equation we have three parameters which are 
fitted. A scale factor on the phase space overlap of the di-pions. The direct production 
D and a. The value for a form the fit is 0.47, which implies r = 0.73 fm. 

We see that the mass shift is a rescattering effect coming from a tttt source which 
is not the p. This tttt source becomes a p° through reinteraction. The mass shift of the 
p seen in jSU], where the p° is photoproduced, is due to the same rescattering effect, 
except that the direct production of the p° in this case is coherent with the tttc source. 
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STAR measured the p° at midrapidity and NA27 measured at the forward region. 
In addition, there is the energy difference of y/s NN = 200 and 27.5 GeV, respectively. All 
these facts are consistent our model, the difference between the p° mass shift measure 
by STAR and NA27 is due to pion multiplicity. Similar arguments is valid for the p 
measurements at CERN-LEP at y/s NN = 90 GeV indicating that there is also significant 
scattering interactions in e + e~ reactions. 

5. Conclusions 

We discussed that the natural mass of the p meson should be measured in exclusive 
reactions only. We showed that a shift on the position of the p° peak has been measured 
before and that the phase space does not account for the p(770)° mass shift measured at 
RHIC, CERN-LEBC-EHS, and CERN-LEP. In addition, we discussed the phase space 
in p + p collisions and concluded that there are significant scattering interactions in p+p 
reactions. These interactions modify the p° line shape in p+p and e + e~ interactions from 
a Breit-Wigner function. Using the rescattering formalism of [2H] , we can reproduce the 
p° mass shift measured by the STAR collaboration in minimum bias p + p collisions. 
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